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We report the recent 4>2 /a results from Belle with KEKB accelerator and Babar with PEP-II 
accelerator. The analysis of B — > tvtv, B — > pp and B — > pn are included in this report. These 
b — > uud decay modes are related to the CKM angle fa jet and the method of <f>2/ot extraction of 
corresponding decays are also included. After combining all the decay modes, the constraint of fa lot 
is (100.2t&°)°. 
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I. INTRODUCTION 

In the Standard Model (SM), CP violating effects in 
the B meson system can be parameterized in terms of 
three Cabibbo-Kobayashi-Maskawa (CKM) p] phase 
angles fa, fa and fa (which can be also written as (3, a 
and 7, respectively). The angle fa/ a can be extracted 
via the b — > uud. Approaches for fa/a extraction from 
the isospin analysis, B — > ttt: and B — > pp, and the 
time-dependent Dalitz analysis, £> — > /97r are reported 
in this document. 

The experimental measurements are from B- 
factories, Belle detector with KEKB accelerator and 
Babar detector with PEP-II accelerator. These two 
detectors are general purpose detectors with energy- 
asymmetric e + e~ accelerators. The Belle detector 
consists of a silicon vertex detector (SVD), a cen- 
tral drift chamber (CDC), an array of aerial threshold 
Cerenkov counters (ACC), time-of- flight scintillation 
counters (TOF), and an electromagnetic calorimeter 
(ECL) comprised of CsI(Tl) crystals located inside a 
superconducting solenoid coil that provides a 1.5 T 
magnetic field. An iron flux return located outside 
of the coil is instrumented to detect if£ mesons and 
identify muons. The Babar detector contains silicon 
vertex tracker (SVT), drift chamber, electromagnetic 
calorimeter, ring-imaging Cerenkov detector (DIRC) 
and a 1.5 T solenoid superconducting magnet. 



II. B — > 7T7T 

In the B° — > tt + it~ decay, the time dependent rate 
is described by 



V q (At) 
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The notation A is used by Belle collaboration which 
is the same as — C used by Babar collaboration. 

The B° — > 7r + 7r - analysis uses 253 fb _1 data col- 
lected by Belle 0. Signal candidates are recon- 
structed by opposite charged tracks which are iden- 
tified as pions and the pion identification is based 
on the combined information from the ACC and the 




= -1 

— Total 

Kti 

continuum 



0.2 0.4 

AE (GeV) 



0.2 0.4 

AE (GeV) 



FIG. 1: AE distributions in the M oc signal region (5.271 
GeV/c 2 < M bc < 5.287 GeV/c 2 ) for B° -> tt+tt~ candi- 
dates with C1Z > 0.86 for (a) q = +1 and (b) q = —1 from 
Belle. 
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FIG. 2: The Belle B° -> tt+tt" At distributions for the 
candidates with C1Z > 0.86 in the signal region (5.271 
GeV/c 2 < M bc < 5.287 GeV/c 2 and \AE\ < 0.064GeV). 
(a) q = +1 and (b) q = — 1. Raw asymmetry, A cp , in each 
At bin with (c) 0<r<0.5 and (d)0.5<r<1.0. The solid lines 
shows the result of the unbinned maximum likelihood fit. 



CDC dE/dx measurements. B meson candidates 
are selected by using the energy difference AE = 



E* B 
M, 



E£ eam and the beam-energy constrained mass 
where EL^ is the CMS 



J 2 



be — yl^bearnJ W B ) ' w 11L1C ^beam 

beam-energy, and E* B and p* B are the CMS energy and 
momentum of the B candidate, respectively. The fla- 
vor of accompanying B meson is identified from the 
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inclusive properties of particles which are not used for 
B° — > 7r + 7r~ reconstruction. To suppress the contin- 
uum background (e + e~ — > qq; q=u,d,s,c), the likeli- 
hood ratio (CTZ) of the event topology based on signal 
MC and sideband data is used for continuum sup- 
pression. The likelihood ratio is optimized separately 
for each flavor tagging quality region. The tagging 
quality is monitored by variable r. After applying 
all above requirement and vertex reconstruction al- 
gorithm, 2820 signal candidates containing 666 ± 43 
7T+7T- signal events (1486 P>° tags and 1334 W tags) 
are obtained. The CP violation parameters S and 
A are determined from the unbinned maximum likeli- 
hood fit to the proper-time difference At distribution. 
They obtained S = -0.67±0.16(stat)±0.06(syst) and 
A = +0.56 ± 0.12(stat) ± 0.06(syst). The correlation 
between S and A is +0.09. Figure and Figure 
show the result from Belle. 



to identify the signal decays. Here pe is the B mo- 
mentum and (Ei, pi) is the four-momentum of e + e~ 
initial state in lab frame. The background suppres- 
sion is based on the angle of sphericity axes and 
Fisher discriminant T formed from the momentum 
flow relative to 7r + 7r~ thrust axis. Unbinned extended 
maximum- likelihood fit is used to extract CP param- 
eters and the likelihood function includes event yield, 
tagging efficiency, mss, AE 1 , T, 9+ , 0^ and At. The 
fit yields S = -0.30 ± 0.17(stat) ± 0.03(syst) and 
C = -0.09 ± 0.15(stat) ± 0.04(syst) from 467 ± 33 
B a — > 7r + 7r _ events. Figure |21 and Figure 01 show the 
result from Babar. 

Using the model-independent isospin analysis 0,0, 
the range of [19°, 71°] for fa/a is excluded by Belle 
with 94.5% C.L. whereas [29°, 61°] is excluded by 
Babar with 90% C.L. 
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III. B 
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FIG. 3: Babar's (a) tues and (b) AE distributions for 
B a — > 7r + 7r~. Solid curves represent the corresponding 
PDFs used in the fit. 




FIG. 4: Babar's At distributions for B° — > tv + tt . (c) raw 
asymmetry for signal events in each At bin. 

The Babar's B° — » 7r + 7r _ analysis Q is based 
on 211 fb _1 data. Signal candidates are recon- 
structed from opposite charged tracks with associ- 
ated Cerenkov angles (6 C ). The particle identifica- 
tion is primarily accomplished by including C in the 
maximum-likelihood fit. The kinematic variables AE 
and M E s = V( s / 2 + Pi ' Pb) 2 /E? - p B 2 are used 



The Babar's analysis of B° — > p°p° |6J| is performed 
with 211 fb _1 data. The reconstruction is made with 
four charged tracks and the particle identification is 
provided by the combining information from DIRC 
and SVT. The continuum suppression is performed 
based on the angle of thrust axis and neural network 
output variable £. Unbinned maximum likelihood fit 
is used to extract the branching fraction by combining 
mES: AE, tt + it~ invariant mass, p helicity angle and 
flavor tagging. From the fit, they obtain the upper 
limit of 1.1 x 10~ 6 at 90% C.L. The m ES distribution 
is shown in Figure [SJ 
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FIG. 5: Projection of B° — » p°p° m,Es distribution ob- 
tained at Babar. 

The time dependent rate of B Q — > p + p~ can be 
parametrized by the same form as Eq. ^ 

The time-dependent CP analysis from Babar is per- 
formed with 211 fb _1 0- Comparing to B — > 7r + 7r~ 
analysis, the B° — > p + p~ need an angular analysis to 
extract the fraction for longitudinal component (/l)- 
An extended maximum-likelihood fit provides /l = 
0.978 ± 0.014(stat)^Q g2g( s y s t) and time-dependent 
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CP parameters, C L = -0.03 ± 0.18(stat) ± 0.09(syst) 
andS L = -0.33±0.24(stat)+°;°4(syst). FigureEDand 
Figure [7| show the result of Babar's analysis. From 
the isospin analysis, Babar obtains 02 jot between 79° 
and 123° with 90% C.L. which is based on the results 
of B — > pp mentioned above, the branching fraction 
of B° ->^P_ @, @ and the results of B+ -> p+ p° 
analysis |10L fll| . 
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FIG. 6: The distribution for the high purity events for 
variables (a) iues, (b) AE, (c) cosine of p helicity angle, 
and (d) m 7r ± n o. The dotted curves are the sum for all 
backgrounds and the solid lines are the total PDF. 




At(ps) 

FIG. 7: The At distribution for the signal enriched sample, 
(a) is B° tagged and (b) is B tagged events, (c) is the 
At raw asymmetry. 

The time-dependent CP analysis from Belle is 
based on 253 fb _1 0. First, the longitudinal com- 
ponent fraction , f L = 0.94ltg;^(stat) ± 0.030(syst), 
and time-dependent parameters, At, = 0.00 ± 
0.30(stat) ± 0.09(syst) and S L = 0.08 ± 0.41(stat) ± 
0.09(syst), are obtained. Figure [8] Figure Eland Fig- 



ure El show the result of Belle's analysis. Com- 
bining other related results 0, 0, Il3j |. Belle obtain 
59° < 2 (cv) < 115° with 90% C.L. 
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FIG. 8: (a) M bc projection in -O.lOGeV < AE < 0.06GeV 
region, (b) AE projection in Mbc signal region. The 
dashed, dotted, dot-dashed, small solid and large solid 
curves show p + p~ + pivir, qq, b — > c, b — > u and the total, 
respectively. 
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FIG. 9: (a) M 7T ± 7T o projection, (b) sum of two cosine helic- 
ity angle distribution. Both plots are inside Mbc — AE sig- 
nal region and satisfy 0.62GeV/c 2 < M^o < 0.92GeV/c 2 . 
The dashed, dot-dashed, dotted and solid curves repre- 
sent p + p~ , pmv, qq + (b — > c) + (b — > u) and the total, 
respectively. 



IV. B 



pn 



Besides B — ► ttit and B — > pp decays, the quasi-two- 
body analysis for B — » pit decay is another candidate 
for the fo/ct extraction. Belle collaboration per- 
formed the time-dependent analysis on B — > p^n^ 
with 140 fb _1 and obtained time-dependent CP 
parameters, C p7r = 



0.25l^ b 7 (stat)+° ;^(sy S t) 



S. 



0.38 



+0.17/ 



-0.28lg;i 

\+0.02/ 



(stat)+o;Jg(syst), 



-0.04 ( s y st ) 



1.08V 

and 



AC, 



P7T 



AS, 



= -0.16 



^_0.18( s * a t) 

-0.30^.23 (stat) ± 0.09(syst), flavor integrated 
charge asymmetry, A p ^ } 
0.02(syst), and direct 
eters, A 
A-+ = 



+0.09 
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(stat) ± 
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CP 
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The branching fraction of B° 



p°ir° is obtained to 
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FIG. 10: The At distribution inside M(, c — AE signal region 
and satisfy 0.62GeV/C 2 < M^o < 0.92GeV/C 2 . (a) and 
(b) are the B° tagged and B° tagged events, (c) is raw 
asymmetry for the events satisfy 0.5<r<1.0. 

be 3.12;° g(stat)±0.33(syst)+o gg(modcl) x 1(T 6 , and 
the CP asymmetry, A CP = -Q.53±g;f£(stat) 
(syst), are obtained by Belle with 357 fb _1 data [lfj. 

Besides the quasi-two-body analysis, Babar per- 
forms the time-dependent Dalitz analysis from B — > 
pn with 7T + 7r~7r° final states ^(| with 192 fb _1 data. 
It is the first direct measurement of fyija, by assum- 
ing isospin symmetry. The decay amplitude can be 
expressed as 

A 37 r = f + A+ + f_A- + f A°, 

A 3 „ = ,f+A + + f_T + f A°. 

The /+, /_ and fo are functions of n w° invariant 
mass that incorporate the kinematic and dynamical 
properties of the B° decay into pir. The Dalitz plot 
distribution for Monte Carlo is shown in Figure ^] 

The time-dependent rate for B° and B tagged events 
are |A^,.(At)| 2 and \A^ 7r (At)\ 2 , respectively. The rates 
are given by 

p -\At\/r B o _ 

|A 3 ±(At)| 2 = — [L4 37r | 2 + MW| 2 

4r B o 

T(\A 3 k\ 2 - L4 3lr | 2 )cos(Am d At) 
±2Im(A 3n A* 3n )sin{Am d At)}. 
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FIG. 11: The Dalitz plot from B° -> n + n~n° Monte Carlo 
without detector simulation. p + iv~ , p~iv + and p°Tv° are 
generated with equal amplitudes. 



There are 27 coefficients for this time dependent 
Dalitz rate, 9 for exponential, 9 for cosine oscillation 
and 9 for sine oscillation terms. Since the B — > p°ir° 
branching ratio is very small, the B — > p°ir° re- 
lated oscillation parameters are fixed to be zero and 
the effect is taken into account by systematics. Fig- 
ure El shows the projection plots for AE' , m^, neu- 
ral network output for continuum background sup- 
pression, At/er(At), m' and 9'. The AE' is the 
transformed AE to deal with the ir° energy depen- 
dence and the o~(At) is the event-by-event error on 
At. The m' and 9' are transformed Dalitz vari- 
ables. The fit yields the direct CP violation values 
A+- = -0.21 ± O.ll(stat) ± 0.04(syst) and A+- = 

-0.47+°;^(stat) ± 0.04(syst). Figure EH shows the 
confidence level of direct CP violation, fo/ct obtained 
from this analysis is (113ti7(stat) ±6(syst))° and con- 
fidence level plot is shown in Figure IT41 



V. SUMMARY 

From the B — > tttt, B — > pp and B — > pn anal- 
ysis, Belle and Babar obtained <\> 2 /a for each de- 
cay separately. The CKM Fitter Group performs a 
global fit by properly averaging all the results and get 
fc/a = (100.21^8°)° [13. Figure US] shows the result 
of this global fit for fa /a. 
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FIG. 12: Distribution of AE' , thes, neural network out- 
put, At /a {At), m' and 9'. The dark, medium and light 
shaded areas represent the contribution from continuum 
events, the sum of continuum events and the B back- 
ground, and the mis-reconstructed signal events, respec- 
tively. 
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FIG. 13: Confidence level contours for the direct CP viola- 
tion. The shaded areas represent 1<t, 2a and 3cr contours, 
respectively. The and ApJ~ here correspond to the 

Belle's A h and A^ , respectively. 




FIG. 14: Confidence level functions for <f>2/ct. The dashed 
horizontal lines corresponds to la and 2a C.L. 
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FIG. 15: The combined fo/u constraint from CKM fitter 
Group. 
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